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1.0 SUMMARY 

The current approach to the Shuttle IMU alignment is based solely 
on the Apollo Deterministic Method (Reference (A)). This method 
is simple, fast, reliable and provides an accurate estimate for 
the present cluster to mean of 1950 transformation matrix. This 
method must remain in the software as the basic estimator, however, 
if four or more star sightings are available, the application of 
least squares analysis can be utilized. The least squares method 
offers the next level of sophistication to the IMU alignment 
solution. The least squares method studied in this working paper 
shows that a more accurate estimate for the misalignment angles is 
computed, and the IMU drift rates are a free by-product of the 
analysis. Core storage requirements are considerably more; 
estimated 20 to 30 times the core required for the Apollo 
Deterministic Method. After the least squares analysis, the next 
level of solution would be the use of the Kalman filter. Kalman 
filtering is not discussed in this paper, but according to Reference (B), 
an entire orbit of star tracker data (90 minutes) is required before 
an accurate solution is computed. Hence, the least squares method 
offers an intermediate solution utilizing as much data that is 
available without a complete statistical analysis as in Kalman 
filtering. 
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• 2.0 IHTRODUCTIOri 

Due to the lack of compensation for certain IMU drift errors, the 
present cluster system slowly rotates away from the desired cluster 
system for which the stored REFSI1AT is valid. Periodically the 
Shuttle IHU's will need realignment to provide a more accurate 
reference to the mean of 1950 coordinate system. With the COAS and 
the Star Tracker systems, line-of-sight (LOS) vectors can be computed 
to known stars on the celestial sphere. This data can then be used 
to determine the validity of the stored REFSMAT and/or v/hen an IMU 
realignment is necessary. The methods discussed in this paper are; 

(1) The Apollo Deterministic Method, requiring tv/o star sightings, 

(2) A direct matrix solution requiring three sightings, (3) and the 
least squares method requiring a minimum of four sightings. 
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3.0 DISCUSSION 

Data from the COAS or Star Tracker is transformed into the actual 

IMU coordinate system (present cluster system) and used in 

conjunction vn’th the M50 (mean of 1950) unit vector to the knov/n 

star (X), The known star unit vector is read from a table based 

on positive identification of the star in the- COAS or Star Tracker 

field of ‘?iew'. The diagram in Figure (1) relates the coordinate 
systems involved with the I^^ alignment method: 



FIGURE (1) 
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The X vector is generated from the measured LOS vector in the Star 
Tracker or COAS coordinate system i>y, 

X»(STSM)Xg, (1) 

where the (STSM) matrix is the Star Tracker or COAS to stable member 
(present cluster) transformation matrix. This matrix is a function 
of the gimbal angle transformation matrix (G) and hence is a function 
of time when the Shuttle is rotating in inertial space. The (STSM) 
matrix is defined from Figure (1) above as: 

(STSM) . C(N)(G)r (2) 

where (N) is a constant matrix relating the Star Tracker or COAS 
to the navigation base reference frame. 

The M50 star vector X from the star table is transformed into the 
ideal (desired) IMU system by the known REFSMAT (R), (a constant 
matrix) as follows: 

Xq = (R) X. (3) 

The vector in the desired cluster system is transformed into 
the actual (present cluster) IMU system by the drift matrix T; 

X = (T) Xq (4) 

The drift matrix is a function of time since it describes the 
relationship between the present cluster and desired cluster 
systems. These two systems drift apart at some small uncompensated 
drift rate o), where the la value for ta is 0.035*’/hr. or about 
0.972 X 10"^®/sec. for the Shuttle IMU's. 
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By substituting equation (3) into (4) 

X = (T)(R) X (5) 

and defining (R‘) = (T)(R), 

the R' matrix is determined when the Apollo deterministic alignment 
calculation is applied to the problem. The solution for the 
instantaneous drift matrix then becomes; 

T = R' R^. (6) 

The misalignment Euler angles may then be extracted from the T 
matrix and displayed to the crew for a criteria to realign the IMU. 


The numerical problem of onboard Shuttle IMU alignment reduces to 
the solution for the drift matrix (T) as a function of time, the known 
star vectors X^, and computed )T^ vectors in the present cluster 
system. Simplification of the diagram in Figure (1) to the basic 
coordinate system involved results in Figure (2). 


(DESIRED 

CLUSTER) 


(PRESENT 

CLUSTER) 



FIGURE (2) 
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An input data table can be made available to the onboard software 

for processing. The format of the input data table is shown in 
Table I. 


POINT 

NO. 

TIME 

M50 

STAR 

LOS 

VECTOR 

COMPUTED 
STAR LOS 
VECTOR IN 
PRESENT CLUSTER 
SYSTEM 

1 

^1 

^1 


2 

^2 

*2 

^2 

3 

^3 

*3 

X 3 

• 

• 

• 

• 

• 

• 

• • 

• * 

• 

• 

• 

• 


t 


I 

n 

n 

n 

n 


TABLE I 


INPUT DATA TABLE 


The computed star LOS vector in the present cluster system is 
calculated using the following procedure: 

Read vertical (5) and horizontal (n) deflection 

angles from the star tracker hare ware system and compute 

the LOS vector (see Figure (1)) as follows; 

£= (1.0 + tan^ 5 + tan^ 

Tnen the components of the vector are. 


( 7 ) 
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Y - tan g 
" l 


( 8 ) 


Xr * tan n 
'2 1 
. - 1 


And the LOS vector above is computed imnediately from 
equation (1) 

X = (STSM) Xg . 


(9) 


3.1 The Apollo Deterministic Solution for IMU Alignment 

When only two points are available, (Table I, N = 2) the 
Apollo Deterministic Method may be used as an estimator for the R' 
matrix (see Figure (2)) and the misaligninent (drift) matrix T. 

This method provides an accurate estimate, for T and R* at the time 
of the observations, but does not estimate the drift. rates, w, of 
the T matrix. 


The method may be described as follows; 

Assuming the first two points of Table I are available, 
create a right-handed system from vectors X^ and X 2 


W = X^ X ^ 


u = U/lUl 
V = V/|V( 
w = 


^u^ v^ w^\ 


X = BX = 1 Vg W 2 IL 

IU3 Vj Wj/ 


( 10 ) 
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Likewise create the same system using and assuming these 
vectors are Invariant with time; 

u* = UVlU'l . 
v» • t/\T*\ 


u* » 

V' * X Xg 

w ' - f, X w* - tTvltT' 


Then form 


u'i v'^ w'^ 

X » CX * 1 u*2 v*2 w '2 . 

iu'3 V3 w*3 


01 ) 


Now solve for t = C^X and substitute this In equation (10); 
X » [Bc"*^] I 


or 


X = [BC^]^ X = (Cb''^) X . 


( 12 ) 

(13) 


Comparing this equation with equation (5) yields; 

R' = Cb'’’ , . (14) 

If, t-j = t 2 (Table I above) the deterministic method will yield an 
exact solution (assuming no star tracker/COAS measurement errors 
and no gimbal angle errors) for the R‘ matrix (equation 14) and 
for the drift matrix from equation (6). However, if t.j f tg, 
geometric errors result in the Apollo deterministic calculation. 

This error results from the fact that the vector (t^) is not 
updated to t 2 (the time of the second observation). A first order 
approximation for this update Is, 
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(t2) * X^ (ti) + (TW Atg) RX^ (15) 

where Atg = t 2 - t, 
and 



The W matrix is the skew-synwietrix angular velocity matrix. Ttie 
oj. 's are the IMU drift rates about each axis. The geometric error 
in the deterministic calculation results from truncating equation 
(15) to, 

(t2) = (t^) . . (17) 

and omitting the resulting error vector 

C2 = (TW At2) R X^ .. (18) 

This error is then further complicated in the calculations forming 
equation (11) and again in equation (14). The final error will be 
directly related to the relative orientation of the two LOS 
vectors and X 2 and the time interval between star sightings, 

i.e.. At = t 2 - t-|. Table II presents a preliminary study of the 

geometric errors based on the 3o drift rates, i.e., tjj= 0.105®/hr 
per axis, with no errors in the star tracker measurements. The 
norms of the errors in the Euler angles are computed to show the 
magnitude of the total error. The differences in the norms of the 
Euler angles for the M50 to present cluster and the misalignment 


^ * • 
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> 

matrix are due to computer round-off mainly in the M50 to present 

cluster computation. The relative magnitude of the Euler angles 

for the M50 to present cluster trans forma cion is +360*. The 

-3 

single precision accuracy expected would then be +1.0 x 10 for 
six significant figures (assuming round-off and slight matrix 
non-orthogonalities). However, the relative magnitude of the Euler 
angles of the misalignment matrix are < ±i.o and the expected 
accuracy is ahout + 1.0 x 10”®. 
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3.2 An Alternate Solution for IMU Alignment 

In the event that 3 points were available In Table I (M « 3), 
the following method offers a solution with possible extension to 
solve for the drift rates. Write equation (4) for each of the data 
points from Table I, 


*1 • T, *0, 
h ' ^2 \ 


Assuming T-j = Tg = T^* I.e., either the drift rates are small or 
the time between sightings is small, ^ 1*^2 = t.^, a matrix can 
be constructed of the form. 

><2 *3) ■ T C^o, ’‘02 V • 

where each vector Is now a column In a matrix and in matrix form 
Is 


(19) 


( 20 ) 


X* Tx, . (21) 

This equation has an lirmedlate solution. 

T » (?)(?/’ . ' (22) 

The restriction here is that x i^JSt be non-singular. Th matrix x 

0 0 

is incleed non-singular if different stars are selected for each of 
the three sightings in Table I. This method, however, has the same 
types of errors as the Apollo deterministic solution. In the 
Shuttle star tracker system, the time between star acquisitions may 
be 5 to 30 minutes, depending on the stars of opportunity. With 
these At's and a 3a IMU drift rate, the geometric error using equation 
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(23) is in the range of 0.01 to 0.1 degrees in the misalignment 
angles.' 

"’’he T matrices of equation (19) are written, 

= T„ + r W At (23) 

noon 

as a first order approximation for the T matrix as a function of time 
and ai. The matrix in equation (23) is defined at some starting 
time t . The time t. is some short interval before t, and At„ = 
t^ - t^. Substituting equation (23) into (19), and adding one 
one more star sighting, we have, 

X, t T„ W At,) 

h (To ^ T„ W Atj) ■ . (24) 

*3 " (To * T„ W Atj) 

*4 = (To ^ T, U At,) . 

By grouping into matrix form, a number of matrix equations can 
be formed. However, a solution for either the T^ matrix or the u 
matrix is not immediately available. This is because the w matrix 
is singular and this first order approximation is valid only for 
small At's. It might be possible to form a closed form approximation 
for T and W using a group of such equations, but further analysis 
of the problem using these methods would be required. The absence 
of a closed form solution for T^ and W suggests the application of 
statistical modeling techniques such as least squares analysis or 
Kalman filtering. 
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3 . 3 Ap plication of Least Squares Analysis for IMU Alignment 

The 1 ine-of-sight vector froip the present cluster as shown in 
Figure (2) and equation (4) for the Nth s'ar sighting is; 

^ <®o. t„) *0 • <“) 

n 

where T„ (9^, u, t^) denotes the T„ matrix which is a function of 
n ' 0 n n 

the misalignment angles 9^, at some time t^, and w represents 
the IMU drift rates. The star tracker measurements can be used to 

A 

compute arid by using an initial estimate for 0^ (from the Apollb . 

A 

deterministic method) and o>, the T^ matrix can be integrated 
(from time t^ to time t^) to obtain an estimate for the vector, 

\ \ <8o* “■ '=n> \ ■ (“> 

n 

By using the difference between the computed measurement Xn and estimated 

vectors X^, . ' 

- X„ = X„ - T„ (e„, w, t„) X„ , ■ (27) 

nnnnn'o n'o^ '' 

n 

a correction to 8^^ and oj can be made by using least squares 
analysis. That is, a correction vector of the form. 



can be computed as a function of the AX^^ vector and applied to an 

iterative solution for 6^ and u; at some time t^. 

0 0 

/A6 \ 

To determine the correction vector i^j . the following analysis 
provides the basic least squares estimator. The vectors X^ are 
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differential of is formed by. 


dx -(5i ro] 

“\36^ 3(0 / Vow / 


(29) 


or 


«n = '’C°) • 


(30) 


The matrix P is the partial derivative matrix whose size is 3n x 6 
(n is the number of star sightings). For small changes, equation 
(29) becomes. 


\3e 3(0 / \A(0 / 


(31) 


The partial derivatives are formed numerically by. 


= ^n - ’^n 

'®o W — 


5(0 


X ft) - X 
^n ^^^5(0 ^ 

5(0 


(32) 


and 


Xn = T (t„) X 

0 on 




( 33 ) 


The estimates for X„ (t_), X (t );-a and X„ (t„). require 

n ' n n ' n'50Q n ' n'5(o ^ 

numerical integration of T (0^, (o) from t^ to t^ the time of star 



WPNo.; E914-8A/B-005 
Page 15 


sighting. The following integrals must be performed: 

T (t^) » / " T (6 q. 0)) dt 
^0 

° ^0 
t 

^ ^®o» “ 

% 

A detailed explanation of the evaluation of these integrals is given 
in paragraph 3.4 of this memo. is then computed by, 

= \ M - T (t„) (35) 

n 

/A0 \ 

and equation (31) is solved for (^1 by 

(to ) , ' 

and 

(to) “ 

This provides the correction vector for equation (28) and with an 
iteration procedure offers a solution to the alignment problem. 

Some preliminary comparison simulations were made with the Apollo 
deterministic method using a maximum of 6-star sighting for the 
least squares solution. Table III shows the comoarison results with 
no star tracker errors, i.e. errors in Apollo deterministic solution 
are those given in equation 18. Table IV shows a comparison with 
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a 100 and 200 arc seconds of random error In the horizontal and 
vertical deflection angles. 


3.4 Integration cf the Transformation Matrix 

In the solution of the IMU alignment method, it is necessary 

to integrate the error matrix relating the actual IMU to the Ideal 

IMU. Since this matrix is the state of the system, state updates 

are required in the analysis. Hence, the matrix must be inte- 

^0 

grated forward ir time to T^. Constant IMU drift rates, are 
assumed. The quantity (t - t^) may take on large or small values 
and the drift rates, are estimated to be about 0.035 degrees 
per hour. The T^ matrix has the derivative; 



(38) 


where the W matrix is the skew-synmetric angular velocity matrix 
(the coordinate axes drift rates, u^.). There are nine elements of 
this matrix each of which must be integrated. Integration of the 
matrix using this form is not practical since truncation errors 
in the numerical process will cause the matrix to become non- 
orthogonal . Integration of the Euler angles is, however possible 
but this method has singularities in the transformation from 
coordinate axis rotation rates to Euler rates. The following 
method utilizes the quaternion rates which are functions of the 
quaternions of the matrix, W, and the coordinate axis rotation rates. 
The quaternion rates are given by 
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= 1/2 (-u^ Qg - ^2 Q3 - U3 94) 

1/2 { 0)^ q, - a>2 q4 + ‘‘*3 q3) 

q3 = 1/^- ( “] q4 + - ^2 q2) 

q4 = 1/2 (-(u^ q3 + 0)2 q2 + <*>3 q^) , 


(39) 


and are non-singular. Using a Taylor series expansion q. 

t+At 

is given by, 

^ ^ ^ ^ + ■ • • + ^ it" . (40) 


Computation of higher order deviatives can be done by using 
equation (39); they are as follows: 

..12 

«li = -4 ^i TT . 


(41) 


where 


*Vf^ = a)-| + ( 1)2 + W 3 

(42) 

II 

1 

1 

(43) 

(4) 14 

‘>1 ■ 16 

(44) 

1 6 

“ii ' 16 "i- 

(45) 

1 -i |6 

<)l=-64‘'l^ 

(46) 

' 6 K« 
Oi ' -64 '’i 

(47) 


and so on for higher order derivatives. By substitution of these 
terms into equation (40) (for n = 7), and simplifying, one obtains 
the following: 



J i' ' 
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= q,* [1.0 - At^ + ^JrTg-y At- - T^JiTg^yAt®] (48) 
t+At t 


Using equation (48) and selecting At such that (uiAt) < 25 (estimated) 
the integral of the transformation matrix can be accurately computed. 
Orthogonality is guaranteed since the matrix will be constructed 
from the final quaternion q. . 

Vat 
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4.0 CONCLUSIONS 

The Apollo deterministic solution provides a basic method for the 
IMU alignment and produces good results even when large star tracker 
errors (random +200 arc seconds in both horizontal and vertical 
deflection angles) are introduced. This method should remain as the 
main computational scheme for IMU alignment. In the event that more 
storage is available during the alignment mission phase, a second 
method such ar the n point least square method (where n i 4) should be 
considered using the Apollo deterministic method output as a starting 
solution. 




WPNo,: E914-8A/B-005 

Page 20 


5.0 REFERENCES 

(A) IMU SOP Section 4. 6. 4. 7. SD76-SH-0013, 4 June 1976. 

(B) H. McQuat, "Further Data on Simulation of a Kalman Alignment 

Filter," CSDL Memo No. lOJ-255-76, 30 June 1976. 


TABLE II 

GEOMETRIC ERRORS IN DETERMINISTIC METHOD 








TABLE III 

COMPARISON OF APOLLO DETERMINISTIC METHOD AND 4 to 6 STAR LEAST 
SQUARES METHOD WITH ND TRACKER ERRORS 
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ERRORS IN MISALIGNMENT ANGLES ARE DUE TO GEOMETRIC ERRORS ONLY 
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COMPARISON OF LEAST SQUARES AND APOLLO DETERMINISTIC METHODS 
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LEO 720 

CEP 1440 6 0.03229 0.06469 .0671 0.1009 



















